Reducing the critical switching current in nanoscale spin valves 
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Abstract 

The current induced magnetization reversal in nanoscale spin valves is a potential alternative to 
magnetic field switching in magnetic memories. We show that the critical switching current can be 
decreased by an order of magnitude by strategically distributing the resistances in the magnetically 
active region of the spin valve. In addition, we simulate full switching curves and predict a new 
precessional state. 

PACS numbers: 73.63.-b, 75.47.-m, 75.70.Ak, 85.75.-d 



The prediction that a spin-polarized current can excite and reverse a magnetization^ 2 - 
has been amply confirmed by recent experiments.—^ The current-induced magnetization dy- 
namics is interesting as an efficient mechanism to write information into magnetic random 
access memories as well as to generate microwaves.— Unfortunately, the critical currents 
for magnetization reversal are still unattractively high.— In this Letter, we apply a previ- 
ously developed microscopic formalism 7 to understand the critical current in spin valves 
quantitatively and propose a strategy to reduce it by up to an order of magnitude. We also 
solve the micromagnetic equations with accurate angle- dependent magnetization torque and 
spin-pumping^ terms and predict switching to a novel precessional state. 

We will first consider a generic F(erromagnetic)|N(ormal)|F spin valve biased by a voltage 
difference V. The two ferromagnetic reservoirs are assumed to be monodomain; the mag- 
netizations differ by an angle 9. Charge and spin currents excited by an applied bias can 
be calculated accurately by magneto-electronic circuit theory 7 - with parameters determined 
from first-principle calculations- that agree well with experimental data.— To this end we 
dissect the pillar into three nodes (the reservoirs and the normal metal) connected by two, 
not necessarily identical resistive elements Gl and Gr. Each of them is characterized by 
the conductances g = + g^, the polarization p = (g** — g^)/g and the normalized mix- 
ing conductance rj = 2g*yg. g'' anc [ gil ar6) respectively, the conductances for electrons 
with spin parallel and anti-parallel to the magnetization and gnJ. is the material parameter 
that governs the magnetization torque. The magnetically active region includes layers of 
thickness up to the spin-flip diffusion lengths from the interfaces. Any resistance outside 
this region is parasitic and not considered here. The conductances are effective parameters 
determined by the resistance of the ferromagnetic and normal metal bulk, that of the in- 
terfaces to the normal metal and the resistance of an eventual outer normal metal that fits 
into the magnetically active region. For simplicity we disregard the bulk resistance of the 
normal metal island and the imaginary part of g^i (fo r metallic interfaces smaller than 10 % 
of the real part2*±i). 

The transverse component of the spin current is absorbed in the ferromagnet 11 and the 
associated spin-transfer torque can excite the magnetization.^ Circuit theory has been used 
to derive analytic expressions for the torques in a symmetrical spin valve as function of the 
angle 9 between the magnetization directions.— However, the spin-transfer torque depends 
strongly on the resistance distribution, even its sign may change with asymmetr y 12 ' 1 ? In this 
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Letter we propose to engineer the spin valve resistance distribution in order to minimize the 
critical current for magnetization reversal from the parallel to the antiparallel configuration 
(the opposite process is much less sensitive to this asymmetry). Expanding the spin-transfer 
torque from circuit theory to first order in 6, the normalized torque i s = \m x (J s x m)\/I c 
on the left magnetization reads 



In an asymmetric structure, Pr — Pl ^ 0, the slope of the spin torque can be considerably 
enhanced compared to the torque in a symmetric structure i s = (h/2e)(p/2)9. 

We investigate a realistic (but non-unique) model for the asymmetry by extending the 
layer sequence from F|N|F to Ni|F|N|F|N 2 . We take here the four N|F interfaces to be equal 
and assume that the ferromagnets are thin enough that the bulk contribution is negligibly 
small. Numerical results require values for the interface resistances that have been measured 
accurately in the current-perpendicular-to-plane geometry . 14 ^ We adopt here Co/Cu inter- 
faces with cross-section 1.26 x 10 4 nm 2 , whence 1/G = 0.0183 Q, p = 0.75 and r\ = 0.38j*^£ 
The asymmetry is modeled by the normal metal sandwich outside the symmetric F|N|F 
structure. The conductivity of the resistive element connecting the left (right) reservoir 
to the adjacent normal metal layer is G\ (G^). The asymmetry is expressed by varying 
the values for G\ and G2 for constant series resistance 1/Gi + I/G2 — 0.37 Q. The total 
collinear pillar resistances can now be calculated by the two current model to be 0.505 and 
0.534 f2, which are typical values for recently fabricated nanopillars.— We assume that the 
right magnetic layer is magnetically hard and is treated as static "polarizer". 

The effective conductance parameters, (g^, Pl and t/l) for the left hand side of the pillar 
consisting of the ferromagnet and the outer normal metal can be calculated in terms of the 
normal metal conductance g\ as 




(1) 



9rVr(Pr ~ Pl) + £lPr(1 - pj) + 9rPl(1 - Pr) 
g L (l ~ Pl) + 9r(1 ~ P 2 r) 



0. 





(2) 



Pl = 



glW - 9 li ) 



2(gi + g n )(gi + g a )9L 
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gi should be replaced by gi to obtain expressions for the right hand side. We parametrize 
Eq. Q as i s = (h/2e)k6, where numerical results for the torque parameter k are given 
in Table for different distributions of the resistance over G\ and (?2- The dependence of 
the torque on G\ : G2 can be understood simply in terms of the spin accumulation in the 
parallel configuration. In contrast to a symmetric structure, it does not vanish in asymmetric 
valves and can have either sign.— The additional spin accumulation that is excited when 9 
becomes finite increases an antiparallel accumulation, and thus the torque, when the extra 
resistance is on the left side. In the opposite case, the initially positive spin accumulation 
is cancelled at a certain angle at which the angular magnetoresistance is minimal and the 
torque vanishes.— 

For our specific example the spin-transfer torque on the left magnetization is enhanced 
by a factor of five when all normal resistance resides on the side of the left magnetic layer. 
The torque can thus be maximized by placing a material with a small spin flip length {e.g. 
platinum) adjacent to the right magnetic layer, as well as a material with a large spin flip 
length {e.g. copper) to the left layer. The magnitude of the torque is enhanced as well when 
all resistance is placed on the other side, but its sign is changed. As shown below, in this 
configuration a reversed current induces switching to a finite angle. 

The magnetization dynamics is described by an extended form of the Landau-Lifshitz- 
Gilbert equation . 17 ' 18 As in previous simulation o 19 ' 20 we adopt a single-domain model, but we 
take into account accurate angle- dependent magnetization torques,— as well as the "dynamic 
stiffness"^ due to spin pumping.— We take the layers to be in the y — z plane and the x- 
axis in the current direction. A uniaxial effective field, -B e ff> and the fixed magnetization 
are chosen parallel to the z-axis. Disregarding dipole and exchange coupling between the 
magnetic layers, both magnetizations in the ground state point along the external field. 
Analytic estimates of the critical current are obtained here by focusing on the instability 
point, at which the current- induced torque exactly equals the damping torque D{6). In the 
presence of in-plane fields, the critical current for complete switching does not necessarily 
agree with the instability point.— We disregard this complication as well as temperature 
induced fluctuations of the magnetizations since they do not interfere with the effect of the 
distributed resistance. 

The maximal viscous damping reads to lowest order in small angles in 9 from the parallel 
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configuration: 

D(9)\ e ^ = aM 1 \B eS \9, (3) 

where a is the Gilbert damping parameter and M\ is the magnitude of the left magnetic 
moment. We obtain numerical results for M\ with the sample cross-section defined above, a 
thickness of 3 nm and a saturation magnetization M s = 1.19 • 10 6 A m -1 . The critical current 
I cc is then given by 



D{9) 



h 



t,L i B (9) 

The total Gilbert damping parameter a consists of cto = 0.006, the bulk Gilbert damping 
parameter and Aa(9) originating from the dynamic stiffness.— ^ We find for Aa| fe0 in the 
limit that the spin currents are efficiently dissipated 



7^ ( tgig n 

8vrMi \ gi + 2g^ 



*«U = ^[^JT + 9 U ). (5) 



The first term in parentheses is the conductance for a transverse spin current from the (left) 
ferromagnet to the left reservoir. The transverse spins escaping to the right are dissipated in 
the ferromagnet when 9 ~ 0, the conductance for these spins is thus rjg/2. The distribution 
of the resistance over G\ and G2 is thus of importance as well for the magnitude of the 
Gilbert damping. Decreasing gi decreases the damping parameter and thus the critical 
current. In Table U the excess damping Aa is given for several resistance distributions. 

The critical currents J CjC can now be calculated assuming an effective field of \B e ff | = 0.2 T. 
We observe that moving the resistance to the side of the switching layer decreases the critical 
current in two ways, by decreasing the excess damping and increasing the torque. For our 
specific model structure the critical current is more than five times smaller in the most 
asymmetric compared to the symmetric pillar. When all resistance resides on the right 
hand side of the pillar, the lowest critical current is achieved by an opposite bias. Not 
only the torque, as shown above, but also the damping is then increased. Measured critical 
currents can be modeled generally well with our model (within 10 %) when anisotropy fields 
are included and, in some samples, the dynamics of the polarizer. 

A Pt layer insertion (with very short spin-flip diffusion length) close to the switching layer 
as fabricated by Kiselev et al., 4 reduces the magnetically active resistance on the left side. 
However, this paper does not report an inverse switching as predicted here. In fact, our 
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calculated critical currents agree best with the experimental ones for a symmetric structure 
without any resistance outside the magnetic layers of equal thickness. This ambiguity might 
be caused by the limitation of our one-dimensional model to accurately describe the three- 
dimensional magnetic polarizing contact in Kiselev et al.s device. 

Finally, we present numerical simulations of complete switching curves. A small initial 
torque is created by starting at 8q = 0.001. In Figure ^ we present the switching curves 
m z (t) of the left magnetization for three resistance distributions. m z (t) is the ^-component 
of the unit vector M\/\M\\. All curves are calculated with c = I c /I c , c — 1-5. We observe 
that for oo : 1 the magnetization switches to an angle between and 7r when the current 
bias is opposite. The origin of this state clearly differs from previously reported precessional 
states,^*^ which required that the applied field is not parallel to the polarizing (fixed) mag- 
netization. It is a direct consequence of the sign change in the torque as function of the 
angle. 12 The switching curves for small deviations from 1 of the different configurations can 
be approximated by (1 — m z {t))/{l — m z (0)) = exp(2«7|.B e fj|(c — l)t). A smaller damping 
parameter thus increases the switching time but decreases the critical current. The angular 
dependence of the spin torque affects the whole switching curve; the torque resembles a sine 
function for 1 : oo, whereas the symmetrical case is closer to Slonczewski's expression.— 

Based on analytic expressions for the spin torque and spin pumping near 9 = in magnetic 
multilayers we conclude that the critical current for magnetization reversal in nano-scale 
spin valves can be reduced by up to an order of magnitude by engineering the resistance 
distribution in the magnetically active region. The spin torque changes sign for specific 
asymmetries giving rise to a new precessional state. After submission of this article, Jiang et 
ali^ reported a strongly reduced switching current by modifying the resistance distribution 
in the nanopillar by a Ru insertion. 

We would like to thank A.D. Kent for fruitful discussions. This work was financially 
supported by FOM, the Norwegian Research Council, and the NEDO joint research program 
"Nano-Scale Magneto-electronics" . 
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FIGURE CAPTIONS 

Figure 1 m z as function of time for switching in magnetic multilayers with different resistance 
distributions. The legend denotes the ratio of G\ : G 2 . 



TABLE CAPTIONS 

TABLE I. The slope of the spin torque at 9 m 0, the increase of the damping by spin pumping 
and the critical current for different asymmetrical configurations of the investigated 
finite element system. 
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FIG. 1: Manschot 



TABLE I: 



Gi : G2 


k 


Aa 


-^c,c 


1 : 00 


0.565 


0.0054 


0.55 mA 


1 : 1 


0.117 


0.0062 


2.86 mA 


00 : 1 


-0.331 


0.0131 


-1.59 mA 



